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ABSTRACT 

X-ray observations of galaxy cluster merger shocks can be used to constrain nonthermal processes in the 
intracluster medium (ICM). The presence of nonthermal pressure components in the ICM, as well as the shock 
acceleration of particles and their escape, all affect shock jump conditions in distinct ways. Therefore, these 
processes can be constrained using X-ray surface brightness and temperature maps of merger shock fronts. 
Here we use these observations to place constraints on particle acceleration efficiency in intermediate Mach 
number (M. ~ 2 - 3) shocks and explore the potential to constrain the contribution of nonthermal components 
(e.g., cosmic rays, magnetic field, and turbulence) to ICM pressure in cluster outskirts. We model the hydro- 
dynamic jump conditions in merger shocks discovered in the galaxy clusters A520 (A4 ~ 2) and IE 0657-56 
(A4 s» 3) using a multifluid model comprised of a thermal plasma, a nonthermal plasma, and a magnetic field. 
Based on the published X-ray spectroscopic data alone, we find that the fractional contribution of cosmic rays 
accelerated in these shocks is < 10% of the shock downstream pressure. Current observations do not constrain 
the fractional contribution of nonthermal components to the pressure of the undisturbed shock upstream. Future 
X-ray observations, however, have the potential to either detect particle acceleration in these shocks through its 
effect on the shock dynamics, or to place a lower limit on the nonthermal pressure contributions in the undis- 
turbed ICM. We briefly discuss implications for models of particle acceleration in collisionless shocks and the 
estimates of galaxy cluster masses derived from X-ray and Sunyaev-Zel'dovich effect observations. 

Subject headings: cosmic rays — galaxies: clusters: general — galaxies: clusters: individual (A520, IE 0657- 
06) — intergalactic medium — shock waves — turbulence — X-rays: galaxies: clusters 
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1. INTRODUCTION 

Astrophysical collisionless shocks are the likely sourc es of 
the observed extra-solar high ene rgy cosmic rays (e.g., iBelll 
ll978t [Blandford & Ostriker n"l978h . Continuous observational 
effort has not yet yielded direct evidence for acceleration of 
cosmic ray nuclei in collisionless shocks, although recently, 
tentative indirect evidence for such acceleration was iden- 
tified in the morph ology of the Tycho supernova remnant 
dWarren et al.l2005l) . and in the high-energy gamma-ray emis- 
sion near the RX J17 13.7-3946 remnant (Ah aronian et al.l 
20071). The theory of d iffusive shock acceleration (e.g., 



Blandford & Eichlen I1987L and references therein) predicts 



that the acceleration efficiency and the spectrum of accel- 
erated particles depend on the Mach number and other pa- 
rameters of the shock (e.g., Giacalone et al. 1997, and refer- 
ences therein). The presence of fossil cosmic rays in the pre- 
shock medium, e.g., from pr evious shocks, can also affect the 
accel eration efficiency (e.g., iKang et aTJl2007t iKang & Jonesl 
120071 and refer ences therein). We here argue that merging 
galaxy clusters are laboratories in which theories of cosmic 
ray acceleration in intermediate Mach number shocks can be 
tested. Observational constraints on particle acceleration in 
such shocks are especially interesting as numerical simula- 
tions suggest that these are the source of a large fraction of 
the co smic rays accelerated in galaxy clusters (e.g. jRyu et alJ 
120031) . 

Fossil cosmic rays, 1 magnetic field, and turbulence may all 

1 Theoretical Astrophysics, Mail Code 130-33, California Institute of 
Technology, 1200 East California Boulevard, Pasadena, CA 91125. 

2 Department of Astronomy, University of Texas, 1 University Station 
C1400, Austin, TX 78712. 

1 Potential acceleration sites of fossil cosmic rays are past accretion and 
merger shocks, giant radio sources, supernovae, and turbulence in the ICM 



contribute to the pressure of the intracluster medium (ICM), 
thereby modifying its hydrodynamic behavior. Such contri- 
bution may alter the interpretation of observations that ig- 
nore them. For example X-ray observations have recently 
been used to estimate shock velocities in two merging clusters 
(see below) neglecting nonthermal components; improved es- 
timates of shock velocities in these and other clusters may 
have to take into account the cosmic ray and other nonther- 
mal contributions to the ICM pressure. Nonthermal pressure 
is also a source of systematic bias when cluster masses are 
estimated from X-ray and Sunyaev-Zel'dovich effect (SZE) 
measurements that are made assuming an hydrostatic equilib- 
rium between gravitati onal forces and thermal pressure gra- 
dients in the ICM (e.g jEnfflin elldll 1 997t iRasia. et al.ll2006t 
Nagai et al. 2007, and references therein). These nonthermal 
biases limit the effectiveness of upcoming cluster surveys in 
the quest to place constraints on the expansion history of the 
universe. 

Evidence for the nonthermal activity in clusters is grow- 
ing. Observed radio and hard X-ray emissions in clusters sug- 
gest a presence of relativistic electrons with Lorentz factor 
of ~ 10 4 . This also suggests a presence of relativistic pro- 
tons that could have been accelerated by the same mechanism 
that has accelerated the electrons. Direct evidence for cosmic 
ray ions in the ICM is, however, still lacking. The nondetec- 
tion in EGRET data of gamma-ray emission expected from 
neutral pion decay in cosmic ray collisions in the ICM (e.g., 
iReimer et al.ll2003l) has so far placed upper limits on the frac- 
tion of cosmi c ray pressure to 5, 20% — 30% in several nearb y 
rich clusters (lEnBlin et al.lll997t iPfrommer & Enfilidl2"004T) . 
Most cluster atmospheres are also substantially magnetized, 
with typical field strengths of order a few /iG out to Mpc radii 

(e.g.,|Berezinsky et all 19971) 
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( Caril li & Tavloril2002l; iGovoni & Ferettill2004l iGovoni et alj 
120061 and references therein). There is likely to be consid- 
erable variation in field strengths (~ 0.1 -40/iG) and topolo- 
gies within clusters. Thus while magnetic fields are likely 
to provide a significant contribution to the p ressure in some 
region s (e.g., along some cold fronts; see IVikhlinin et al.1 
2001), it is yet unclear what is the average energetics in 
the cluster outskirts. Numerical simulations of cluster for- 
mation also suggest that subsonic gas motions (turbulence) 



Norman & Bryan 1999; Ricker & Sarazin 


2001; Nagaietal. 


20031: iFaltenbache 
20061 iNaeai et al.1 


r et al.ll2005UDolae et alJl2005l iRasia et all 
20071). Further investigations of nonther- 



mal phenomena in clusters are hence critical for the success 
of upcoming X-ray and SZE cluster surveys, as our ability to 
estimate cluster masses hinges on a precise characterization 
of the nonthermal components. 

In this work, we show that shock waves that form dur- 
ing merging of galaxy clusters can provide unique con- 
straints on nonthermal processes in clusters. Recent Chan- 
dra X-ray observations have revealed that shock waves with 
Mach numbers M k 2-3 that form during the merging 
of galaxy clusters are accompanied by distinct X-ray sur- 
face brightness and temperature discontinuities. To date, 
a density and temperature jump along axis of symmetry 
of a cluster merger bow shock has b een recovered in the 
cluster A520 (iMarkevitch et al l 120051) an d in IE 0657-56 
dMarkevitch et al J2002t lMarkevitchl2006l) . The relative mag- 
nitudes of the density and temperature jump depend on the 
contribution of nonthermal components (e.g., cosmic rays, 
magnetic field, turbulence) to the pressure of the undisturbed, 
upstream ICM. They also depend on the efficiency of particle 
acceleration in the shock, on the escape of the accelerated par- 
ticles from the shock, and on the amplification of turbulence 
in the shock. 

Here we make the first attempt to use the X-ray observa- 
tions of merging clusters to constrain their particle accelera- 
tion and the fractional contribution of nonthermal components 
to the pressure budget of the ICM. In §|2] we model the hydro- 
dynamic jump conditions and Mach number using the multi- 
fluid approximation. We model the effect of cosmic rays, a 
tangled magnetic field, and turbulence, on observed gas jump 
conditions in the shock. In § [3] we utilize this model to de- 
rive constraints on nonthermal components in merger shocks 
in galaxy clusters A520 and IE 0657-56. We show that cur- 
rent observations can limit the efficiency of particle accelera- 
tion, and that the fractional contribution of nonthermal fluids 
to ICM pressure may be constrained with future, improved, 
X-ray spectroscopy and SZE observations. In §|4] we discuss 
implications of these results in context of particle acceleration 
models and for the estimates of galaxy clusters masses with 
X-ray and SZE observations. In §[3] we summarize our main 
conclusions. 

2. LIMITS ON PARTICLE ACCELERATION IN MERGER SHOCKS 
AND NONTHERMAL PRESSURE IN ICM 

2.1. Shock Jump Conditions 

The ICM fluid consists of a thermal component (electrons 
and ions in mutual thermal equilibrium) and number of non- 
thermal components (nonthermal cosmic rays and magnetic 
fields comoving with the thermal component). Turbulent gas 
motions and electromagnetic waves can also contribute pres- 
sure to the ICM. We ignore the electromagnetic waves in the 
current treatment. We also temporarily ignore turbulence, but 



return to discuss its role on qualitative level in § 12.81 The 
ICM fluid can be modeled with an adiabatic equation of state 
with an effective adiabatic index 7. 2 In what follows, we use 
indices "u" and "d" to denote the shock upstream and down- 
stream, respectively. Assuming cylindrical symmetry, conser- 
vation of density, momentum, and energy fluxes across the 
shock dictates jump conditions on the symmetry axis of the 
bow shock that read 

/O u V u = PdVd, 



^u + PuVu = -Pd + PdV'd, 



2 7u 



7u 



1 



Pu)=Vd 
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7d 



1 



(1) 



where p, is the mass density of the thermal gas (the mass den- 
sity of the nonthermal particles is negligible), v, is the fluid 
velocity, and P, is the total pressure (z = {u, d}) 

The pressure is a sum of electronic, ionic, cosmic ray, and 
magnetic field contributions, Pj = P e l +Pi,, - + Pcr,/ + Pb,i- The 
effective adiabatic index 7,- equals 
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where ccr,; = PcR,i/Pi an d eg.; = PB,i/Pi are, respectively, the 
fractional contribution of cosmic rays and magnetic fields to 
the total pressure, while e n u = ecr,; + cb,/ is the total fractional 
contribution of nonthermal pressure components. The tem- 
perature jump across the shock of the thermal electrons is then 



7e.d 

T 

1 en 



(l-e„t,d)[27u/(7u-D- 



•1] 



(l-£nt.u)[27d/(7d-D- 
where r = pa/p u is the compression ratio. 



•1] 



(3) 



2.2. Thermal Electron Relativistic Corrections 

Thermal ions are nonrelativistic in cluster shocks and thus 
7y = 5/3. The adiabatic index of thermal electrons may 
differ from this value because of relativistic corrections. 
The adiabatic index of thermal electrons equals the ratio of 
the rest-frame pressure to internal energy density, j e = 1 + 
Pe/£e, where £ e is the internal electron energy density (e.g., 
AchterbergHaDEHl 



1 J^47T P \p 2 /c 2 + mlr l / 2 Up)dp 
7e + 3 f °° 4tt P 2 [(p 2 c 2 + m 2 e c 4 y/ 2 



where 



f e (p) oc exp 



m e c 2 ]f e (p)dp 

(p 2 c 2 + mlc A ) l l 2 ' 
k%. 



(4) 



(5) 



is the thermal electron momentum distribution. Defining 
the electron temperature in units of the electron rest energy 
6 e = kT e /m e c 2 , the adiabatic index in equation (01 can con- 
veniently be expressed i n terms of modifie d Bessel functions 
of the second kind (e.g., Kunik et al. 20031) 



7e,j = 1 + e ,/ 



3e e ,+ 



gi(e-j) 
W-j.)' 



(6) 



2 We define the adiabatic indices via 7 = 1 +P/e, where e is the inter- 
nal energy density. The adiabatic index thus calculated may differ from 
ainP/ainp. 
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In the limit 9 £ < 1, relevant to cluster merger shocks, the adi- 
abatic index is 7 e = 5/3-(5/6)e e + £>(9 2 ). Note that when 
kT e i or kT & u becomes comparable to m e c 2 , the right side of 
equation (Q depends on the electron temperatures. For ex- 
ample, if the nonthermal pressures are ignored (i.e., e nt ,■ = 0), 
a relativistic correction of about 10% applies to the expected 
temperature jump and inferred Mach number for conditions 
similar to those in the bow shock in IE 0657-56 (T u = 10 keV 
and r = 3), while a correction of only about 1 % applies for the 
lower-temperature and weaker shock in A520 (T u = 5 keV and 
r = 2.3). 

2.3. Fossil Cosmic Rays 

The undisturbed ICM may contain an intracluster pop- 
ulation of fossil cosmic rays that could have been pro- 
duced in the high Mach number accretion shocks (Ai ~ 
10-100), in previous merger shocks, in active galactic 
nuclei, and in starburst-associated phenomena (see, e.g., 
Berezinskv et al.lll997t iFuiita & Sarazirj|l200U iMiniati et ail 
2001b iGabici & Blasill2003b ISarazinll2004l) . After a merger 
bow shock passes a fluid element in the ICM, the fluid el- 
ement will contain the original fossil cosmic rays, some of 
which may have been further accelerated in the shock. The 
shock may also accelerate new cosmic rays drawn from the 
thermal pool and accelerated in the shock for the first time. 

The cosmic ray adiabatic index depends on whether their 
pressure is dominated by Newtonian or relativistic particles, 
which in t urn depends on the deta ils of the cosmic ray spec- 
trum (e.g., Achterbe rg et al.| [l984). The effective cosmic ray 
adiabatic index lies in the range 4/3 < 7cr < 5/3; the elec- 
tron cosmic rays are mildly or fully relativistic while the pro- 
tons that dominate the cosmic rays pressure can be Newtonian 
or relativistic. Thus, assuming that electron and proton cos- 
mic rays are close to mutual equipartition, we expect 7cr to 
range between 7cr ~ 4/3 (for relativistic electrons and pro- 
tons) and 7c r ~ 13/9 (re lativistic electrons and Newtonian 
protons; e.g.,|Kpnigl 1980). While in principle 7cr may dif- 
fer between the upstream and the downstream, when there is 
a significant population of upstream cosmic rays, it is reason- 
able to expect that these will also dominate the downstream, 
and thus that 7cr is approximately constant across the shock. 

In addition to being reaccelerated by shock acceleration 
mechanism, fossil cosmic rays may be accelerated adiabat- 
ically during the shock compression of the magnetic field. 
However, shock compression may very well be nonadiabatic 
as the Larmor radius of a mildly relativistic electron in a mi- 
crogauss field, tl ~ 10 y cm, may be much larger than the tran- 
sition layer over which the fluid density jumps. 3 Therefore, it 
is possible that cosmic rays do not gain energy during com- 
pression. Nevertheless, the compression increases the cosmic 
ray number density and therefore its pressure by at least a fac- 
tor of r. Thus we express the fractional cosmic ray pressure 
in the downstream as 



ecR.d = r [ — | £cr,u + £acc 



(7) 



where the first term is the pressure of the fossil cosmic rays 
if the energy of individual cosmic rays remains unchanged 
as they pass through the shock, and the second term is any 

3 If pressure across the transition layer is not mediated by cosmic rays but 
by plasma instabilities, the width of the transition layer will be of the order 
of the plasma skin depth of the thermal protons, which is S ~ 10 cm. On the 
other hand, adiabatic acceleration will take place when rL <S 5. 



additional pressure of cosmic rays that are accelerated in the 
shock. The acceleration term includes pressure due to accel- 
erated particles from the upstream thermal pool that are added 
to the nonthermal pool, and also due to fossil nonthermal par- 
ticles that are reaccelerated by a shock acceleration process 
and/or by adiabatic compression. 

2.4. The Cluster Magnetic Field 

The effect of the magnetic field on the shock depends on 
its strength and topology, both of which are constrained very 
poorly, especially in the cluster outskirts. Faraday rotation 
measurements suggest that the magnetic pressure can reach 
~ 10% of the thermal pressure in some clusters or regions 
within clusters, and that in others it is less than 1% (e.g., 
Carilli & Taylor! 120021: iGovoni & Ferettil 12004 iGovoni et al.1 
2006). 4 Given the large uncertainty in the field strength, we 
treat eg as a free parameter in what follows. 

We assume that the magnetic field of the undisturbed ICM 
is tangled and isotropic on scales relevant to the hydrodynam- 
ics of the cluster merger shock. This approximation allows us 
to use the unmagnetized form of the jump conditions, equa- 
tion (Q]i, to describe the effect of magnetic field on shock hy- 
drodynamics. The approximate coherence length of the field 
in some galaxy clusters, ~ 10 kpc, is shorter than the curva- 
ture scale of the merger shock, ~ 100 kpc, so that one can 
average over the fluctuating field orientation near the shock. 5 
Furthermore, the measurements of density and temperature 
jumps across the shock are based on deprojection of the X-ray 
map assuming a cylindrical symmetry of a bow shock — thus 
an averaging of observables on the shock curvature scale is 
implicit in the reported shock jump data. 

Even if the upstream field is isotropic, the downstream field 
may be anisotropic as a result of a preferential amplification 
of the magnetic field component that is perpendicular to the 
shock, B±. In this case the perpendicular component will still 
be isotropic in the plane of the shock, and can thus be param- 
eterized by the ratio of the average parallel and perpendicular 
field energy densities b = (Bj_)/2(B|). The effective adia- 
batic index of the magnetic fluid, 7b, which is defined as the 
parameter that correctly quantifies the behavior of the mag- 
netic pressure in equation and may depend on the orien- 
tation of the magnetic field, can be expressed in terms of the 
averages of the components 7^™ of the electromagnetic part 
of the energy-momentum tensor in the shock frame via 



1b 



i7> 



(8) 



Here, (3 is the shock velocity in units of the speed of light, and 
r = (1 -/3 2 )" 1 / 2 is the Lorentz factor of the shock. Compo- 
nents of the energy-momentum tensor read 



1 
1 



T™ = ^- [(2r 2 - 1 )B\ - fijj ] , 



With these, for a nonrelativistic shock we obtain 



7b : 



4b 
2b+l' 



(9) 



(10) 



4 A magnetic field that is tangled on very small scales is not well con- 
strained by the Faraday rotation measurements. 

5 The same approximation was employed by Markevitch et al. ( 2005). 
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The value of jb in equation ( fTOb equals 7# = 2 in a perpendic- 
ular shock (b = 00) and vanishes in a parallel shock (b = 0). 
If the field is isotropic on average in the shock upstream 
(b = 1), we recover the value jg.u =4/3. For an isotropic up- 
stream field, the downstream field will have b > 1 and thus 
4/3 < 7B . d < 2. 

Relating es,d to □ requires an understanding of the shock 
structure. But even if the shock structure is unknown, we can 
limit £ba m two extreme cases, when field generation in the 
shock itself (as may take place in unmagnetized shocks) is 
neglected. Assuming that magnetohydrodynamic jump con- 
ditions apply within the shock transition, and that the par- 
allel and the perpendicular fields do not transform into each 
other, then Bu j = Bii u and Bx,d = rB± u , implying that e^.d = 
(2r 2 - l)(P u /Pd)£B,u- Th e other extreme assump tion, which 
was previously made by Mar kevitch et alJ (120051) . is that the 
field is isotropic and remains isotropic throughout the shock 
compression, in which case tB.d = r 4 ^(P u /Pd)eB,u, and thus 
the magnetic field behaves as a relativistic or photon gas that 
is adiabatically compressed in the shock. We expect eg.d to 
lie between these two limits, assuming that no new field is 
generated within the shock. 

2.5. Escape of the Highest Energy Cosmic Rays 

Another process that may affect the shock and that we have 
ignored so far is the escape of the highest-energy cosmic rays 
that are accelerated in the shock. This process may take place 
if the shock efficiently produces a hard cosmic ray spectrum, 
and the highest-energy cosmic rays can es cape and remove 
energ y from the shock transition (see, e.g.. [Achterber get alJ 
1984). This energy leakage can be parameterized by the frac- 
tion Q of the incoming energy flux that escapes the shock. 
The leakage affects the shock jump conservation relations 
such that the left-hand side of the energy flux conservation 
relation — the last of equations (Q]) — must be multiplied by the 
factor 1 — Q. With this, equation ([3]) generalizes into 

= (1 - £nt,d)[2(l - 07u/(7u - 1) - r-' - 1 + Qr/(r- 1)] 
(l-£„t. u )[27d/(7d-l)-r-l + e>-7('--l)] ' 

(11) 

The effect of escape is similar to the effect of particle accel- 
eration, as it removes a fraction of the incoming energy from 
the downstream thermal fluid, thereby reducing the tempera- 
ture jump. Note, however, that if the measured value of r is 
below the one expected in a thermal shock, the deviation from 
the thermal temperature jump can only partially be attributed 
to the escape, since when e acc = 0, the value of Q must by 
definition be zero. 

2.6. Sensitivity of the Temperature Jump 

We have one relation (equation [3] or its generalized form, 
equation [TTTi and several unknown parameters describing the 
upstream and downstream contribution of the cosmic rays and 
magnetic field to the fluid pressure. Without making assump- 
tions about the nature of the nonthermal fluids, given a mea- 
surement of the density jump r and the temperature jump r, 
we can place constraints in the joint parameter space spanned 
by these parameters, but cannot recover the fractional pres- 
sures themselves. Additionally, when only soft X-ray data 
are available (as from Chandra and XMM-Newton), the un- 
certainties in t typically exceed the uncertainties in r. There- 
fore, before we proceed to explore the joint parameter space, 
we discuss the effect of the variation of various parameters on 



r, assuming that r is accurately measured and is held fixed. 
The temperature jump is most sensitive to the fractional pres- 
sure in cosmic rays accelerated or reaccelerated in the shock, 
e acc . A nonzero value of e acc reduces the value of r, e.g., in 
the ideal case in which the upstream nonthermal pressure van- 
ishes, e nt , u = 0, we find that r drops by about a factor of two 
between e acc = and e acc = 0.3 (r = 2-3). The reason for this 
sensitive dependence is that the production of accelerated par- 
ticles in the shock saps a fraction of the incoming energy flux 
out of the thermal component of the downstream, thereby re- 
ducing r. Increasing Q results in a similar effect, for similar 
reason, on r. 

The effect of changing 6cr, u and eg.u on the temperature 
jump is more subtle since a high nonthermal pressure in the 
shock upstream implies that there will also be a high nonther- 
mal pressure in the downstream; it is the balance between the 
two that determines r. Therefore, r depends only weakly on 
the upstream nonthermal components. Increasing ecr.u while 
keeping the rest of the parameters constant always increases r 
slowly. For example assuming that the magnetic pressure van- 
ishes (e^.u = 0) and that 7cr = 4/3, r increases by a factor of 
« 1 .2 between 6cr. u = and £cr, u = 0.3 for r = 2.3 and e acc = 0. 
Having a constant e acc = 0.3 reduces the change in r to a factor 
of w 1 .1. Increasing cb , u also mostly results in a slightly larger 
temperature jump. If the perpendicular field is strongly am- 
plified in the shock, 7^ jumps in the shock, thereby increasing 
the effective downstream adiabatic index 7,5, and thus increas- 
ing r. If on the other hand we treat the field as a relativistic 
gas, then eg , u and r are positively correlated for e acc = and 
are weakly anticorrelated for e acc 0. 15. 

In conclusion, r is strongly anticorrelated with e acc and 
is typically weakly positively correlated with £cr,u and e^.u- 
Therefore an accurate measurement of t and r can tightly con- 
strain e acc . If the measured value of r falls below that expected 
in a purely thermal hydrodynamic shock given a precise mea- 
surement of r, then a nontrivial lower and upper bound can 
be placed on e acc , but such a constraint cannot be obtained for 
£cr,u and e^.u- If, on the other hand, the measured value of r 
is higher than that expected in a purely thermal hydrodynamic 
shock, only an upper limit on e acc can be placed, i.e., e acc = 
remains viable. Then, however, the measurement places a 
lower limit on the upstream nonthermal pressure. However, 
based on the measurement of r and r alone, one cannot sep- 
arate the partial contributions of the cosmic ray and magnetic 
components. 

2.7. Implications for Shock Velocity Estimation 

The velocity of cluster merger shock is of great interest 
since it is a stepping stone toward relating the dynamics of 
the ICM to th e dynamics of the dark matter in galaxy cluster 
mergers (e. g., Hava shi & Whi te 2006; Farrar & Rosenll2006l : 
iMilosavrievic et al.ll2007hlSprin gel & Farrar||2007|) ."Typically 
the shock Mach number is inferred from the compression ratio 
r, which can be measured relatively accurately in X-ray maps, 
under the strict assumption that all pressures are thermal (the 
observed constraints on r are typically used for consistency 
check with this assumption). Taking the nonthermal pres sure 
into account (but ignoring cosmic ray escape; see § 12.51 , the 
Mach number is given by 

M 2 = 2r 7u/(7u-l)-27d/(7d~l) {m 

7g,u( 1 - e„tu)( 1 - r- 1 ) [2 7d /( 7d - 1 ) - r- 1 ] ' 

where 7 g u = 5/3 is the adiabatic index of the upstream ther- 
mal gas, for which relativistic corrections are negligible. This 



NAKAR ET AL. 



5 



relation implies that high nonthermal pressure in the upstream 
(downstream) increases (decreases) M.. For example, assum- 
ing that e nt u = and cosmic rays are rather efficiently acceler- 
ated in the shock, e acc = 0.15, the inferred value of the Mach 
number must be revised by » 10% -20% downward of the 
value inferred for a thermal shock, when r is in the range 2-3. 

If, in addition to r, the value of r is accurately measured, the 
constraints that can be placed on the Mach number are tighter, 
since the upstream and downstream nonthermal components 
are no longer entirely free. In this case, when r and r are 
held constant, a larger nonthermal pressure results in a higher 
Mach number. For example if the nonthermal component be- 
haves as relativistic gas (7„ t = 4/3) and r and r are related as 
they would be in a purely thermal shock, t = (4r— l)/r(4 — r), 
the shock Mach number equals 



M 



: _ 3(4r - 1)(1 - e n t,u)/(4 - r) + 5r+ 3re nt , u - 8 
7u(l-£nt,u)(l-r- 1 )(7-r) 



(13) 



in which case A4 increases by « 10% for e nt . u = 0.3, compared 
to the purely thermal shock, for r in the range 2-3. 

2.8. Turbulence 

In addition to the cosmic rays and the magnetic field, tur- 
bulence also contributes pressure to the ICM. Turbulence in 
the ICM is expected to be driven by gravitational clustering 
(accretion and merging), and by outflows associated with ac- 
tive galactic nuclei. Hydrodynamical simulations of galaxy 
cluster formation in which the ICM is treated as an ideal fluid 
universally demonstrate that turbulent pressure in the ICM is 
non-negligible, and that its contribution to the total pressure is 
an increasing function of radius from the center of the cluster. 
The fractional turbul ent pressure measured i n the simulations 
is e tur ~ 0.06-0.36 dNorman & Brvanll999b. e t „ r ~ 0.05-0.1 
dRicker & Sarazinl2001l) e, m ~ 0.04-0.09 dNagai et al.l2003l: 
Faltenbacher et a l. 2005), and e tur ~ 0.05-0.3 dDolag et al l 



2005). 6 Spatially-resolved gas pressure maps obtained from 



XMM-Newton observations of the Coma galaxy cluster show 
a scale-invariant pressure fluctuation spectrum on the scales 
of 40 to 90 kpc, which was analyzed to place a lower limit on 
the fr actional turbulent pressure of e tur > 0. 1 (Schuecker et al. 
l200l . 7 

The physics of the interaction of a shock wave with a tur- 
bulent upstream is complex and poorly understood even in 
unmagnetized, ideal fluids. The primary theoretical uncer- 
tainties are the amount of amplification of turbulence in the 
shock and the affect of turbulence on the shock structure. The 
shock trans ition becomes nonplanar in the presence of turbu- 
lence (e.g., Rotman 1991), and this nonplanarity affects the lo- 
cal shock jump conditions that short-wave length fluctuation s 
experience while crossing the shock (e.g.. IZank et al.l 2002). 
Different analytical approximations (e.g., the rapid distortion 
theory and the local interaction analysis) and direct numer- 
ical simulations do not a lways agree with each other (e.g., 
lAndreopoulos et ail l2000, and references therein). Therefore 
we do not attempt to include turbulent pressure in our quanti- 
tative calculations. However, since turbulent pressure can af- 
fect the observed jump conditions, we discuss it qualitatively 

6 Hydrody namic simulations of intracluster turbulence w ere also recently 
carried out by Fujita et al. 12004). Subramanian et al. (2006), and Vazza et al. 
(MS) . 

7 Turbulence in the ICM can be detecte d in other clusters given an 
X-ray detector with high spectral resolution tlnogamov & Sunyaev] 120031 : 
Sunyaev et al. 2003). 



using a simple model (lLeld ll992) based on the rapid dis- 
tortion theory applied to homogene ous turbulence (see, e.g., 
Batc helorf[T953HJacquin etaill 19931) . 

iLeld (1 1992b derives the averaged density, momentum and 
energy conservation equations of ideal fluid in the shock 
frame (his equations 9-11). Assuming homogeneous tur- 
bulence and a cylindrically symmetric distribution of turbu- 
lent fluctuations, the conservation equations are reduced to the 
form of equations (HJ with a turbulent pressure and effective 
adiabatic index of 



7tur = 



3 + 2fr,- 
l+2b,' 



(14) 



Here, the notation is such that any fluctuating quantity / is 
decomposed in two ways / = /+/' = /+/", with / denoting 
the average value of the quantity, f=pf/p denoting the mass- 
weighted average, and /' and /" denoting the corresponding 
fluctuating parts, while as before, i = (u,d). Just as we did for 
the magnetic field (see § 12.41 ). we use cylindrical symmetry to 
parameterize the turbulent field with a single parameter = 

v" -Vj_ ; /(2vj' j .vj' j ), which in the isotropic case equals unity. 

The temperature jump r and the shock Mach number M. 
can be calculated for a given upstream turbulent pressure frac- 
tion etur,u = Ptx,\x/Pu an d anisotropy parameter b u if the ampli- 
fication of the turbulence in the shock is known. 8 ILeld (1 1992b 
uses the rapid distortion theory to derive the shock amplifica- 
tion. This approximation assumes that the mean turbulent flux 
amplitudes are much smaller than their mean flow counter- 
parts, that turbulent fluctuations cross the shock much faster 
than the corresponding eddy turnaround times, and that the 
mean flow does not vary much on the length scale of an eddy. 
In particular, the otherwise planar shock transition is assumed 
to have not been distorted, and rendered nonplanar, by the 
fluctuations. In this theory, the par allel and pe rpendicular tur- 
bulence is amplified according to (lLeld ll992) 



! II u 



l 



— r 
4 



6r-l 



1 (r 2 -l) 3 / 2 



tan 



Vr 2 -1 



1- 



r+1 



1 



1 (r 2 - 1) 3 /2 



tan 



(15) 



which in the case of isotropic upstream turbulence, b u = 1 and 
7mr,u = 5/3, yields effective adiabatic index and downstream 
pressure fraction for turbulence 

7tur,d ~ TT~ H i (I 6 ) 



llr+4' 



Ctur.d ! 



etur,u(6r- 1) 



5t(1 -e tur „) + e tU r,u(6r- 1) 



In this simple model of amplification of turbulence in the 
shock, the effect of a pre-existing turbulent component is sim- 
ilar to that of the pre-existing cosmic rays and magnetic pres- 
sure, namely, the upstream turbulence is only weakly, posi- 
tively correlated with r and A4. For example, assuming that 

8 In case of isotropic turbulence e lur = ygML^/Q+ygM^), where jVftur 
is the turbulence Mach number and 7 g is the thermal gas adiabatic index. 
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turbulence is the only nonthermal component, the tempera- 
ture jump increases by a factor ps 1.15 and Ai increases by 
a factor « 1.3 between e tur u = and e tur ,u = 0.3 for compres- 
sion ratio r = 2.3. Very similar results are obtained if instead 
of the rapid distortion theory, we employ the line ar interac- 
tion analysi s to calculate the jump conditions (e.g.. lLee et al.1 
fl99l fl997h . Therefore, we conclude that for a weak, sub- 
sonic, unmagnetized turbulence of the ICM, the effect of tur- 
bulence on the thermal gas temperature jump of the shock is 
similar to that expected in the presence of other nonthermal 
components. 

3. RESULTS 
3.1. A520 

iMarkevitch et al.l (|2005) analyze a 67 kilosecond observa- 
tion with Chandra ACIS-I of the bow shock in the galaxy 
cluster merger A520 at z=0.203 and estimate the density and 
temperature jump along the axis of symmetry of the shock. 
The upstream and downstream temperatures are, respectively, 
T u = 4.8!j;| keV and T A = 1 1 keV, while the density jump 
is r = 2.3 ± 0.3 (all errors are at the 90% confidence levels). 
These values are consistent with a shock Mach number of 
A4 « 2. We use equation (0 to explore the constraints that 
can be placed on the acceleration of particles in the shock 
and on the fractional pressure in the nonthermal components 
in the shock upstream, and equation ( fTTT i to constrain energy 
leakage from the shock in A520. 

The predicted downstream temperature in the case of a 
pure thermal gas, = 10 keV for T u = 4.8 keV and r = 2.3, 
is consistent with the measured value. However, given that 
the expected temperature is below the median measured tem- 
perature, little room is left for significant particle accelera- 
tion in the shock. To place constraints on e acc , we first take 
CB i = 2 = 0, and carry a Monte Carlo search in the remain- 
ing parameter space (ecR,u,£acc)- We draw a large set (10 5 ) 
of the observed parameters (density and temperature jump) 
from the observed distributions. 9 For each set of observed 
values we find all the combinations of (ecR, u ,eacc) in the do- 
main (0 < ecR.u < 0.3,0 < e acc < 0.25) that are compatible 
with the observations. For each point in the (ecR, U) face) plane, 
we calculate the number of instances that the corresponding 
shock is compatible with the generated "observed" values of 
r and r. The resulting number, properly normalized, provides 
the Bayesian likelihood of the shock being characterized by a 
given pair (ecR,u, £acc), assuming a uniform prior in the domain 
considered here. 

Figure [TJ shows the probability distribution for (ecR.u,£acc) 
for the relativistic cosmic rays (7cr = 4/3). The two contours 
divide the plane so that the cumulative probability constrained 
above each is 0.33 and 0.05 of the total. The figure shows that 
e acc < 0.1 at 95% confidence levels for any value of £cr, u , and 
that current observations do not provide a constraint on £cr,u- 
Carrying out a similar analysis for relativistic electrons and 
Newtonian proton cosmic rays (7cr = 13/9), which would be 
expected if the cosmic ray pressure were dominated by parti- 
cles with typical energies of 10 1 - 10 4 times the thermal en- 
ergy, yields an upper limit of e acc < 0.15. Finally, the limit on 
the efficiency of the acceleration of Newtonian cosmic rays, 
7cr = 5/3, as expected if the cosmic rays are accelerated only 
to several times the thermal temperature, is e acc < 0.25. Natu- 

9 We approximate each observed distribution by two half-gaussians that 
peak at t he median observed valu e and satisfy the 90% confidence range re- 
ported by Markevitch et al. (20051). 
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FIG. 1. — The probability distribution of (ecR,u> e acc) for relativistic cos- 
mic rays in the bow shock in the galaxy cluster A520, assuming £.b,w 

= and 

e CR,u < 0.3. The contours divide the plane so that the cumulative distribution 
above the contour includes only 0.33 (lower contour) and 0.05 (upper con- 
tour) of the total probability. The probability is calculated via a Monte Carlo 
simulation (see text). 
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FIG. 2. — The marginalized probability distribution of e a cc for several 
different nonthermal contributions to the gas pressure. The probability is 
calculated via a Monte Carlo simulation (see text) for various ranges of flat 
priors on ecR.u an d e fl.u, and various levels of magnetic field amplification 
in the shocks, as indicated in the legend. The cosmic rays are relativistic 
(7CR =4/3), unless noted otherwise. The solid and dashed lines contain 90% 
of the distributions (the rest is in dotted line tail). This figure shows that the 
limits on e acc are insensitive to the assumptions and the priors that we choose 
for the upstream nonthermal components. 



rally, allowing for Q > when deriving the constraints on e acc 
yields a tighter upper limit on e acc . The minor effect of vary- 
ing upstream magnetic field and upstream cosmic ray pressure 
on the limits that can be placed on e acc is explored in Figure 

m 

Repeating the same analysis for various values of < e# u < 
0.3 while assuming 6cr, u = 0, or taking a constant 6cr, u and 
a flat prior on 6b, u in the same range, for the two limiting 
forms magnetic field behavior at the shock (B^a = rB±_ u and 
; see § |2.4t , yields a similar upper limit on e acc . 
The reason for the weak impact of magnetic pressure on the 
constraints that can be placed on particle acceleration is the 
strong dependence of r on e acc and its weak dependence on 
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£b,u- Therefore, just as we found for £cr,u, current observa- 
tions do not provide a significant constraint on ee, u . 

We also carried out the same analysis assuming no particle 
acceleration, e acc = 0, but allowing for an escape of cosmic 
rays from the shock, in order to constrain energy leakage from 
the shock, Q. Although this scenario is artificial (since for 
face = we also expect Q = 0), this analysis provides an upper 
limit to the value of Q. We find that the observation of A520 
limit Q < 0.1 in its merger shock. 

The dependence of the temperature jump in equation (01 on 
ccr.u and e£ U is weak. Therefore current observations do not 
put significant constraint the presence of a relativistic com- 
ponent in the pre-merger ICM. Figure [TJ shows that a future, 
improved X-ray spectroscopy across the bow shock, such as 
with a longer exposure with Chandra or XMM-Newton, may 
exclude the purely thermal scenario, e acc = ecr.u = 0. If mea- 
sured value of r falls below the thermal prediction, our anal- 
ysis will yield a lower limit on e acc to accompany the current 
upper limit. If the measured value of r falls above the ther- 
mal prediction, the measurement will imply a lower limit on 
ecR.u- For example, if we artificially reduce the present mea- 
surement uncertainties in T u , T<i, and r by a factor of 3 while 
assuming that the mean values of these observables remained 
unchanged, the data would require e acc < 0.05 and e nt u > 0.05 
at 95% confidence levels. However decoupling e ntu into its 
cosmic ray and magnetic field components cannot be accom- 
plished given a measurement of r and r alone. 

From the den sity jump data al o ne an d assuming a purely 
thermal shock, Ma rkevitch et alJ (120051) estimate the Mach 
number of t he m erger shock in A520 to be Ai =2.1^ 3. As we 
discuss in § 12.71 the likely presence of a non-negligible non- 
thermal pressure requires a modification of the Mach num- 
ber estimate as in equation (fT2l : adding a nonthermal pres- 
sure in the upstream increase the Mach number. Figure [3] 
shows the Mach number probability for several scenarios in 
which a nonthermal fluid is present with a fractional contri- 
bution to the pressure that is allowed by the present observa- 
tions. The true value of the Mach number depends on the 
fractional nonthermal pressure and the amplification of the 
magnetic field and turbulence in the shock. For example, for 
ccr,u = £b,u = 0.15 and assuming significant magnetic field am- 
plification in the shock, the true Mach number can be as high 
as M « 2.7 

3.2. IE 0657-56 

Gas properties across the merger shock in IE 657-56 (the 
"bulle t" cluster at z=0.296) were measured by Mar kevitchl 
(2006) using a 500 kilosecond observation with Chandra 
ACIS-I. They find a density jump of r » 3, which corresponds 
to a Mach number of M. ~ 3. The measured temperatures 
are T a k, 9 keV and a lower limit > 32 keV at la confi- 
dence levels. The shock in IE 0657-56 is stronger than that 
in A520 and is thus more propitious for detecting particle ac- 
celeration. Unfortunately, the high downstream temperature 
Td complicates an accurate measurement of the temperature 
jump with the high resolution X-ray telescopes Chandra and 
X MM-Newto n 

Mar kevitchl ((2006) does not report the errors on some of 
the measurements and therefore we cannot quantitatively con- 
strain the presence of nonthermal components in the shock 
upstream and downstream. Figure H] shows the model predic- 
tion for T<± as a function of 6cr. u and e acc for the case of rela- 
tivistic cosmic rays, 7cr = 4/3, assuming that eg.u = 2 = and 
taking r = 3 and T u = 9 keV. The figure shows that constraints 
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FIG. 3. — The Mach number probability distribution based on the observed 
density and temperature jumps across the shock in A520 for several different 
contributions of nonthermal components to the gas pressure. The solid lines 
contain 90% of the distributions (the rest is in dotted line tails). In all cases 
the cosmic rays are relativistic and the upstream magnetic field is isotropic 
(7CR = 7s.u = 4/3). The legend indicates the fractional upstream nonthermal 
pressure and the shock amplification of the field in each case. 

that can be made in IE 0657-56 are qualitatively similar to, 
although less stringent than, those in A520. The minimum 
value of the temperature jump r allowed by the measurement 
is high and barely consistent with a purely thermal shock. It 
does not leave much room for particle acceleration; we tenta- 
tively infer e acc < 0.15. 

The nonthermal components may affect also the Mach 
number of the bullet cluster merger shock, which recently 
stirred a discussion a bout its compatibility w i th standard cos- 
molo g ical models ( Hayashi & White 2006; Farrar & Rosenl 
2006; Milosavl jevic et al.1 120071; ISpringel & Farrarl |2007|) . 
Markevitchl d2006l) finds M = 3.0 ±0.4 , which corresponds to 
r = 3.0!» 23, assuming no nonthermal contribution and neglect- 
ing relativistic corrections due to the high electron tempera- 
ture. Allowing for cosmic rays with pressure up to equiparti- 
tion, < ccr.u < 0.3 and eb = 0, while requiring temperature 
and density jump consistent with observations, > 20 keV 
and 2.77 < r < 3.17, yields the limits 2.3 < M < 3.7. The 
lower bound is obtained for r = 2.71, e acc = 0.07, ecr.u = 0, 
and 7d = 20 keV. The upper bound corresponds to r = 3.17, 
e acc = 0, e C R.u = 0.3, and T d = 45 keV. 

4. DISCUSSION 

4.1. Particle Acceleration in Collisionless Shocks 

Current measurements place constraints on the efficiency 
of particle acceleration in cluster merger shocks. This of- 
fers a very unique opportunity to test models of particle ac- 
celeration in astrophysical collisionless shocks. Estimating 
the acceleration efficiency in collisionless shocks is a diffi- 
cult problem, which is severely complicated by the fact that 
it remains unknown which among a number of possibilities 
is the primary acceleration mechanism. Even in the lead- 
ing candidate mechanism, the diffusive shock acceleration, 
estimates of acceleration efficiency range widely because of 
a number fundamental theoretical uncertainties, concerning 
the fraction of thermal p articles that are injecte d into the 
acceleration process (e.g., Malkov & Druryll200ll and refer- 
ences therein), the nonlinear influence of the accelerated par- 
ticles on the hydrodynamic profile of the shock wave (e.g., 



8 



SHOCK CONSTRAINTS ON NONTHERMAL PRESSURE IN CLUSTERS 



0.25 



0.2 



0.15 



o 
o 

CO 

CO 



0.1 



0.05 



T =9keV ; r=3 






u 

























135 



H30 



25 

< 



20 



15 



0.05 



0.1 0.15 

Si. 



0.2 



0.25 



0.3 



'CR,U 



FIG. 4. — T& as a function of ecR.u and e a cc (7cr = 4/3) for r = 3 and 
7d = 9 keV, which correspond to the values reported for the bullet cluster IE 
0657-56. The 32 keV and 20 keV cont ours are the la and 2cr lower limits 
on 7j as measured by Markevitch 1 2006) (the 2ct limit is a rough estimate 
derived from the plots in Markevitch 2006). 



Drury & V6Ddll981t lAchterberg et al.lll984t iGiacalone et all 
1997t iKang et al.l l2002l I2007I) . and the ampli fication of the 
magnetic field by plasma instabilities (e.g.. iLucek & Bell 
2000 [ iBell & LuceMl200U ISchlickeiser & Shuklal 120031: iBelll 
2004tl2005HSchekoch ihin et a l.ll2005HMedvedev et alj|2006l 
Therefore, any estimate of particle acceleration efficiency is 
strongly affected by the specific assumptions and approxima- 
tions employed in a self-consistent shock model. Here we 
do not attempt to carry out a detailed comparison of our re- 
sults with various scenarios for particle acceleration. But, to 
demonstrate the power of the constraints that can be obtained 
from merger shock dynamics, we discuss our results in the 
co ntext of the predic tions of the particle accelera tion model 
of IKang et all (120071 see also lKang & Jonesll2007l) . 

A variety of investigations point to a strong dependence of 
acceleration efficiency on the shock Mach numb er, whereby 
strong er shocks produce higher efficiencies. Kan g et all 
(2007) report an investigation of acceleration efficiency be- 
havior in diffusive shock acceleration simulations in quasi- 
parallel shocks with a Bohm diffusion coefficient, a self- 
consistent treatments of particle injection from the thermal 
pool into the accelerat ion process, and A lfven wave propaga- 
tion. In their Figure 5. lKang et al.l (12007) plot the dependence 
of the acceleration efficiency r/(Ai) on the shock Mach num- 
ber M, where 77 is defined as the energy flux in downstream 
cosmic rays, divided by the bulk kinetic energy flux in the up- 
stream medium entering the shock. In the limit e acc -C 1, the 
parameter rj in Kang et al. is related to our e acc via 



3 (M 2 + 3)(5M 2 -l) 



10 



M 4 



(face < 1). 



(17) 



Thu s, we have 77 sa 2 5 e acc for M = 2 and 77 ss 2e acc for M = 
3. 10 |Kang et al. (2007) detect a strong dependence of 77 on the 
presence of preexisting cosmic rays in the shock upstream, 
ccr.u- The strong dependence can be attributed to inefficient 
injection at low Mach numbers in their model. For M. = 2, the 
parameter 7/ jumps from to 0.15 as ecR, u increases from to 
0.23. In stronger shocks, for M. = 3, the parameter 77 jumps 
from 0.1 to 0.25 for the same increase in £cr, u - 

10 We have here ignore d the modification of shock Mach number by the 
nonthermal pressure, see § 12.71 



Converting our current limits on e acc into limits on r\, we 
find 77 < 0.2 for M « 2 and 77 < 0.3 for M » 3. These limits 
are ma rginally consistent with the predictions of Kan g et aT] 
(2007) for any a ssumed value of £cr,u- However, if as 
Kang et al. (2007) argue, e acc is itself a sensitive function of 
ecR.u, then improved measurements of e acc that can be ob- 
tained with additional observations with existing X-ray tele- 
scopes can provide tighter limits on £cr,u- Moreover, a posi- 
tive measurement of e acc can provide a model-dependent con- 
straint on 6cr,u- 

4.2. X-ray and SZE Cluster Mass Estimates 

Galaxy cluster surveys in which the cluster masses are 
measured accurately can be used as powerful cosmological 
probes of dark matter and dark energy. The mass estimates 
are plagued by systematic uncertainties that must be under- 
stood and quantified before the requisite mass measurement 
accuracy is achieved. Nonthermal pressure due to cosmic 
rays, magnetic fields, and turbulence, is a source of a sys- 
tematic bias when cluster masses are estimated on the ba- 
sis of the assumption of hydrostatic equilibrium between 
gravitation al forces and thermal pressure gradients in the 
ICM (e.g.. lOstriker et al.l2 005: Ras ia. et al.l2006tlNagai etail 
2007, and references therein). The hydrostatic mass profile of 
a spherically-symmetric cluster is given by 



M(< )- ^ ( dP s + dP rt 
Gp g \ dr dr 



(18) 



where M{< r) is the mass enclosed within radius r, while P g 
and P nt are the thermal and the nonthermal contributions to 
the pressure. The thermal gas provides a significant fraction 
of the total pressure support, and this pressure is measured 
directly with current X-ray and SZE observations. The con- 
tribution of the nonthermal pressure, on the other hand, is cus- 
tomarily assumed to b e relatively small ( < 10%) outside of a 
cluster core (see e.g.. iNagai et al.l [20071) . and it is often ig- 
nored in the hydrostatic mass estimates based on X-ray and 
SZE data. However, present observations do not yet constrain 
the nonthermal pressure in the regime in which it dramatically 
affects the calibration of the hydrostatic mass estimates. If 
not accounted for, these nonthermal biases limit the effective- 
ness of upcoming X-ray and SZE cluster surveys to accurately 
measure the expansion history of the universe. Detailed inves- 
tigations of the sources of nonthermal pressure in clusters are 
thus critical for understanding their effect on the properties 
of the ICM and the utility of clusters as precision cosmologi- 
cal probes. We proceed to discuss how future observations of 
cluster merger shocks will have the potential to place unique 
constraints on the nonthermal pressure in the unshocked ICM, 
thereby improving cluster mass estimates. 

4.3. Prospects for Future Constraints of Nonthermal 
Pressure 

As discussed in ^3] current measurements alone do not 
place strong constraints on the presence of a nonthermal com- 
ponent in the unshocked ICM in both systems. However, the 
improved constraints that can be obtained with existing X- 
ray telescopes can provide useful lower limits on nonther- 
mal pressure and their effects on the X-ray and SZE clus- 
ter mass estimates. In the case of the shock in the cluster 
merger A520, the current constraints, which are based on a 
67 kilosecond observation with Chandra, can be improved 
significantly with follow-up observations with Chandra or 
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XMM-Newton. Therefore, we identify this system as the most 
promising one in which our method may yield a positive 
shock-hydrodynamic detection of a nonthermal component. 
While the stronger shock in the cluster merger IE 0657-56 
may be more efficient at accelerating particles, it will be more 
difficult to improve the measurement of the shock temperature 
jump in IE 0657-56. This is because the very high tempera- 
ture of the shock downstream medium (Tj « 30-50 keV) lies 
far outside of spectral sensitivity window of X-ray telescopes 
with arcsecond resolution, which in turn renders it difficult to 
measure the temperature jump in the narrow post - shock layer. 
Hard X-ray observations (e.g., iPetrosian et al.l 120061) with 
RXTE, Integral, or Suzaku, combined with a model of ICM 
fluid flow (e.g.. iMilosavlievic et alj|2007t Upringel & FarraH 
l2007h . may help to pin down the downstream temperature. 
Alternative intriguing possibility is that future high resolution 
SZE observations may be able to measure the downstream 
temperature on the basis of the relativistic SZE, which should 
be prominent in the high temperature downstream. 

In the next few years, gamma-ray observations of galaxy 
clusters may provide tight constraints on the fractional con- 
tribution of nonthermal particles to the pressure of the ICM. 
Assuming that gamma-ray emission from the decay of neu- 
tral pions is the primary emission channel, measurements 
with the new gamma-ray telescope GLAST can be used to 
place population-averaged li mits on the hadronic cosmic-ray 
pressure support in c lusters (lAndo & Nagai I I2007L see also 
iBerrington & Permed l200l iBlasi et al.ll2007l) . These forth- 
coming constraints, combined with improved X-ray spec- 
troscopic measurements of the jump conditions in merger 
shocks, may enable a separation of the upstream nonthermal 
pressure into its constituent components. 

Finally, comparisons of the hydrostatic mass estimates with 
those derived from gravitational lensing observations can, in 
principle, provide an important handle on nonthermal biases. 
Note that this approach is not practical for individual clusters, 
because lensing measures the mass in a projected aperture that 
cannot be directly compared to the mass within a sphere of the 
same radius, to which the hydrostatic mass is sensitive. But, 
it might be possible to compare different estimators in an av- 
erage sense, while accounting for the effe cts of asphericity o f 
clusters and projection effects (see § 5.2 in lNagai et al.ll2007h . 

5. CONCLUSIONS 

We model the effect of particle acceleration and nonther- 
mal pressure components on shock jump conditions in non- 
relativistic shocks. We focus on intermediate Mach number 
shocks, with Mach numbers in the range M. = 2- 3. We ap- 
ply this to the merger shocks in galaxy clusters A520 and IE 



0657-56 and place the first constraints on the efficiency of 
particle acceleration in these shocks. Our main results are as 
follows. 

1 . The temperature jump of the thermal gas in the shock 
depends strongly on the efficiency of shock particle accelera- 
tion. Efficient acceleration can reduce the temperature jump 
by more than a factor of two for a constant compression ratio 
in the range r = 2— 3. 

2. The correct effect of nonthermal pressure in the up- 
stream, such as fossil cosmic rays, magnetic field, and tur- 
bulence, on the shock jump observed in the thermal gas can- 
not be derived at this point, because we lack an understanding 
of the interaction of these components with the shock. How- 
ever, for a wide range of reasonable assumptions and analytic 
approximations, we find that nonthermal pressure in the un- 
shocked ICM has only a minor effect on the downstream tem- 
perature (at a fixed compression ratio), and that in general, a 
high upstream nonthermal pressure increases the temperature 
jump in the thermal gas. 

3. The combination of strong dependence of the tempera- 
ture jump on particle acceleration and weak dependence on 
upstream nonthermal pressure enables derivation of mean- 
ingful constraints on the efficiency of particle acceleration in 
cluster merger shocks, even with current observations. Future, 
more accurate X-ray and SZE observations of these shocks 
may yield meaningful constraints on the upstream nonthermal 
pressure as well. 

4. Nonthermal pressure and shock particle acceleration can 
also affect the Mach number that is inferred from the observed 
compression ratio r by tens of percent. When the temperature 
jump is poorly constrained, the Mach number is antic orrelated 
with efficient particle acceleration and is positively correlated 
with upstream nonthermal pressure. 

5. In the two observed high contrast galaxy cluster merger 
shocks, A520 and IE 0657-56, we constrain the efficiency 
of acceleration of relativistic particles to be e acc < 0.1 and 
e acc < 0.15, respectively. We find that considerable upstream 
pressure can increase the Mach number of the shock in A520 
to reach M. ~ 2.7, much higher than the inferred value 2.1 ob- 
tained assuming an absence of nonthermal components. The 
true Mach number of the shock in IE 0657-56 can be in the 
range 2.3 < Ai < 3.7 with the compression ratio of r = 3 al- 
lowing for nonthermal pressure components. 
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